Magnetism in doped two-dimensional honeycomb structures of III-V binary 

compounds. 
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Using first-principles plane-wave calculations systematic study of magnetic properties of doped 
two-dimensional honeycomb structures of III-V binary compounds have been conducted, either for 
magnetic or nonmagnetic dopants. Calculations show, that all cases where magnetic moment is 
non-zero are energetically more favorable. For such cases band structure and (partial) density of 
states were calculated and analyzed in detail. The possible applications of these structures were 
also discussed. 
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I. INTRODUCTION 

Since its discovery in 2004 graphene [l| draws much at- 
tention beacuse of uniqe features of this two-dimensional 
system. Graphene is composed of a sp 2 -bonded carbon 
atoms forming honeycomb structure. It has very interest- 
ing electronic structure with characteristic, linear energy 
dispersion near K point of Brilloin zone. Summary of the 
subject can be found for example in 
Shortly after, experimental techniques allowed fabrica- 
tion of other new two-dimensional materials, like BN 
and M0S2 honeycomb structures Q or ZnO monolay- 
ers @. The discovery of such stable two-dimensional 
material like graphene triggered search for similar struc- 
tures made from different compounds. Up to now many 
of these hypotetical structures constructed from silanene 
(2D Si) and germanene (2D Ge) 0,i, III-V compounds 
1, SiC or ZnO Q have been studied theoretically. 
On the other hand graphene and other nano-scale ma- 
terials are recognized as future building blocks of new 
electronics technologies [l(|, including spintronics (e.g. 
11]). In the case of low (one- and two-) dimensional 
structures problem arises because of famous Mermin- 
Wagner theorem [l2|], which prevents ferro- or antifoer- 
romagnetic order to occur in finite temperatures, which is 
essential for spintronics and other modern applications. 
This started the theoretical and experimental search for 
magnetism in graphene and other structures. One of the 
most promising directions is emergence of magnetism in 
such structures as an effect of presence of local defects 
[l3j . According to works of Palacios et al. [l4[ and, inde- 
pendently, of Yazev [l5[ single-atom defects can induce 
ferromagnetism in graphene based materials. In both 
cases, the magnetic order arises as an effect of presence 
of single-atom defects in combination with a sublattice 
discriminating mechanism, in agreement with Lieb's the- 
orem 17]. Based on these findings several theoretical 
studies have been conducted in search for magnetism in 



low-dimensional structures either for graphene and BN 
18] or other (hypothetical) structures like SiC Q- 
In this paper influence of local defects on magnetic struc- 
ture of two-dimensional honeycomb structures of GaN, 
A1N and InN have been analysed by means of ab-initio 
calculations. Since bulk versions of these compounds 
are very important semiconductors in todays electron- 
ics it would be interesting to check whether such two- 
dimensional materials could have non-zero magnetic mo- 
ment. Despite of the fact that neither of them have been 
yet synthesized, calculated cohesion energies @ suggest 
that such structures would be stable and their experi- 
mental procurement is highly probable. 
Next section contains computational details followed by 
results. Last section concludes this work. 
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II. COMPUTATIONAL DETAILS 



To investigate magnetic properties of GaN, A1N and 
InN honeycomb structures a series of ab-initio calcu- 
lations have been conducted with use of DFT VASP 
code [jj, [U with PAW potentials [Hj]. For 
both spin-unpolarized and spin-polarized cases exchange- 
correlation potential has been approximated by general- 
ized gradient approximation (GGA) using PW91 func- 
tional [22j]. Kinetic energy cutoff of 500 eV for plane- 
wave basis set has been used. Supercells of size 4x4x1 
have been checked to be large enough to prevent de- 
fects interact with its periodic image. In all cases for 
self-consistent structure optimizations, the Brillouin zone 
(BZ) was sampled by 20x20x1 special k points. All struc- 
tures have been optimized for both, spin-unpolarized and 
spin-polarized cases unless Feynman-Hellman forces act- 
ing on each atom become smaller than 10 -4 eV/A. A 
vacuum spacing of 12 A was applied to hinder the inter- 
actions between monolayers in adjacent supercells. Cal- 
culated lattice constants are in agreement with |6j. 
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III. RESULTS 

As mentioned, non-magnetic honeycomb sheets can 
attain spin polarized states due to presence of local 
defects. In this work two kinds of defects have been 
analysed - vacancies and subsitutions. 
For all three compounds a vacancy was generated first 
by removing a single atom, Al, Ga, In or N from each su- 
percell, then the atomic structure was optimized. In all 
cases structures with single N vacancy are non-magnetic, 
while Al, Ga or In vacancies induce non-zero magnetic 
moment, equal to 3.00 (is- This is in disagreement with 
Lieb's theorem, which states that magnetic moment 
should be equal to 1.00 /is when one of sublattices has 
exacly one atom more/less than the other (e.g. Nai - 
= ± 1). This discrepancy can be addressed to charge 
transfer from Al(Ga,In) to N. Fig. [T] shows density of 
states (DoS) for spin-polarized Al-vacant A1N, on which 
difference between majority spin (up) and minority spin 
(dn) in the vicinity of Fermi level (horizontal line) can be 
observed which is the main source of non-zero magnetic 
moment. Analysis of calculated partial magnetisation 
shows that almost all magnetic moment is situated on 
p-states of N atoms located in the area of vacancy. This 
is in full agreement with previous studies of vacancies in 
SiC Q. 
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FIG. 1: Density of states of Al-vacant A1N. 

In the case of substitution the procedure was as follows. 
For all three compounds various single foreign atoms have 
been substituted, then structure has been optimized. In 
the case of A1N and GaN, Al or Ga has been subsituted 
by atoms from 4-th period of periodic table from K to 
Zn, including Na and Mg for A1N. In the case of InN, 
In has been subsituted by atoms from 5-th period of pe- 
riodic table from Rb to Cd, excluding Tc. In all three 
compounds, N has been subsituted by C, B and P atoms. 
Table I shows calculated magnetic moments and differ- 
ences between total energy of spin-unpolarized and spin- 
polarized states, AE = E„ sp - E sp , i.e. positive value of 
AE means that spin-polarized state is more energetically 
favored. This is the case for all compounds with non-zero 
magnetic moment - the largest energy differences (up to 
1.25 eV) are for instances with highest values of magnetic 



moment. For A1N highest values of induced magnetic mo- 
ment are in case of Mn, Co (4.00 /i#) and for Fe (4.26 
/is). Similar situation can be observed for GaN doped 
with Mn, Co and Fe. With decreasing number of d-shell 
electrons value of magnetic moment drops as well as for 
the case of Zn which has d-shell closed. In the case of 
InN this tendency holds although values of magnetic mo- 
ments are much smaller, being the highest for Ru and Rh. 
Figs. [3] and [4] show mechanism of generation of magnetic 
moment in the case of GaN doped with Ni. Left plot 
of Fig. [3] shows bandstructure of Ni-doped GaN vs. un- 
doped one (which is a semiconductor with bandgap equal 
to 2.30 eV, calculated within GGA) in spin-unpolarized 
case. One can see the formation of doping bands in the 
vicinity of Fermi level. The top right plot shows band- 
structure of GaN+Ni in the spin-polarized case, where 
can be observed quite large splitting of these bands be- 
tween spin up and down bands. 

Top plots of Fig. @]show density of states for spin po- 
larized case. Right one shows total DoS vs. DoS of spin 
up and down. Fermi energy is almost exactly in the mid- 
dle of splitted up and down DoS. Left plot shows vicinity 
of Fermi energy more closely, where large splitting of up 
and down DoS can be observed. Since almost all elec- 
trons occupying vicinity of Fermi level are d-shell elec- 
trons, which can be read from partial density of states 
(bottom left plot) the mechanism of magnetic moment 
emergence becomes clear. Calculations show that this 
mechanism is universal for all structures having non-zero 
magnetic moment doped with transition metal elements. 
In case of doping with alkali metal elements and alkaline 
earth metal elements only in Na- and K-doped A1N cal- 
culations show non-zero magnetic moment (1.88 and 1.70 
/US) respectively for Na and K). Mechanism of formation 
of magnetic moment is similar to the case of vacant struc- 
tures (since Na and K have only one valence electron). 
Fig. [5] shows DoS for spin-polarized Na-doped A1N, on 
which difference between spin up and spin down in the 
vicinity of Fermi level can be observed. 
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FIG. 2: Density of states of Na-doped A1N. 

This is similar to situation depicted on Fig. [TJaltough 
in the case of Na-doped structure splitting is smaller. 
In case of substitution of N atom by C, B and P, only C- 
doped structures had non-zero magnetic moment, which 
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was equal to 1.00 /is in all compounds. Magnetic states 
in all cases are lower by about 0.1 eV than nonmagnetic 
states. 

IV. CONCLUSIONS 

Ab-initio calculations have been conducted for vacancy 
and subsitution defects in honeycomb A1N, GaN and InN 
compounds. Calculations show that in all three com- 
pounds vacancy of Al, Ga or In respectively gives mag- 
netic moment of 3.00 ug, which is interesting conclusion 
from application point of view. On the other hand substi- 
tution of Al or Ga by transition metal elements (Mn, Fe, 
Co) can give even higher value of magnetic moment (4.00 



/is). Since techniqe of implantation of metal atoms into 
2D surface has been recently reported [23| , it is also very 
promising direction. On the other hand substitution by 
non-metallic atoms or substitution of nitrogen atoms by 
IV or V group atoms does not give significant magnetic 
moment. Thcs results may give a hint for experimental- 
ists seeking for two-dimensional magnetic materials. 
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FIG. 3: Bandstructure of Ni-doped GaN. Details in text. 
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FIG. 4: Density of states for Ni-doped GaN. Details in text. 
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AIN 

for Al Na Mg K Ca Sc Ti V Cr Mn Fc Co Ni Cu Zn 
h{hb) 1-88 0.00 1.70 0.00 0.00 0.96 0.00 2.98 3.99 4.26 3.96 2.98 1.85 0.00 
AE (eV) 0.06 0.00 0.04 0.00 0.00 0.11 0.00 0.70 0.98 1.09 0.94 0.38 0.06 0.00 



GaN 

for Ga K Ca Sc Ti V Cr Mn Fc Co Ni Cu Zn 

h{Hb) 0.00 0.00 0.00 0.90 1.72 2.99 3.99 4.51 3.97 2.97 1.32 0.00 

AE (eV) 0.00 0.00 0.00 0.12 0.07 0.66 1.15 1.25 0.91 0.30 0.01 0.00 

InN 

for In Rb Sr Y Zr Nb Mo Tc Ru Rh Pd Ag Cd 

h{/j,b) 0.00 0.00 0.00 0.00 0.99 0.39 - 1.38 1.91 0.42 0.00 0.00 

AE (eV) 0.00 0.00 0.00 0.00 0.16 0.01 - 0.04 0.10 0.01 0.00 0.00 



TABLE I: Magnetic moments and total energy differences be- 
tween spin up and spin down states for different transition 
metal elements dopants. 



AIN 

for N B C P 

0.00 1.00 0.00 
AE (eV) 0.00 0.14 0.00 

GaN 

for N B C P 
fj,{n B ) 0.00 1.00 0.00 
AE (eV) 0.00 0.11 0.00 

InN 

for N B C P 
fj,{n B ) 0.00 1.00 0.00 
AE (eV) 0.00 0.10 0.00 



TABLE II: Magnetic moments and total energy differences 
between spin up and spin down states for different dopants. 



